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Abstract 

The electronic structure of the perovskite LaCoOa for different spin states of Co ions was calculated in 
the LDA+U approach. The ground state was found to be a nonmagnetic insulator with Co ions in a low- 
spin state. Somewhat higher in energy we found two intermediate-spin states followed by a high-spin state at 
significantly higher energy. The calculation results show that Co 3d states of t2g symmetry form narrow bands 
which could easily localize whilst eg orbitals, due to their strong hybridization with the oxygen 2p states, form 
a broad a* band. With the increase of temperature which is simulated by the corresponding increase of the 
lattice parameter, the transition from the low- to intermediate-spin states occurs. This intermediate-spin 
(occupation tisej) can develop an orbital ordering which can account for the nonmetallic nature of LaCoOa 
at 90 K<T<500 K. Possible explanations of the magnetic behavior and gradual insulating-metal transition 
are suggested. 

1 INTRODUCTION 

Among the systems showing semiconductor to metal transition LaCoOs is especially interesting due to the 
fact that is also displays very unusual magnetic behavior, often associated with low-spin (LS) - high-spin (HS) 
transition Although a large number of investigations have been carried out since the early 1960's the character 
of the transition and the nature of the temperature dependence of the spin state is still unclear. For example 
the temperature dependence of the magnetic susceptibility shows a strong maximum at around 90 K followed 
by a Curie- Weiss-like decrease at higher temperatures Q| which was interpreted by the authors as a LS to a HS 
transition. The semiconductor to metal-like transition occurs in a range of 400 ~ 600 K, well above this transition. 
The insulating nature of LaCoOs below 400 - 600 K was attributed by Raccah and Goodenough to an ordering 
of LS and HS Co'^^ ions in a NaCl like structure , with the itinerant electrons in a broad band formed by the 
transition metal Cg orbitals. This structural change however has not been observed in crystallographic studies. 
Recent neutron scattering experiments |^ suggest that the semiconductor - metal transition is not predominantly 
magnetic in origin but a clear picture of what its nature is absent. 

Photoemission (XPS) and X-ray absorption (XAS) studies on the other hand have been interpreted in terms 
of a spin change at the semiconductor - metal transition [Q. There are however several discrepancies in the 
interpretation. The Co 2p XAS line shape at room temperature looks like LS Co"^+ which is not consistent with 
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a spin state transition at 90 K. However, the valence band XPS spectrum at 300 K is quite different from that 
of a HS compound. So the interpretation of the two transitions is still uncertain. 

In most of the previous studies one has assumed a rather ionic and ligand field- like starting point. In this 
picture the possibility of intermediate spin (IS) states as well as local symmetry lowering, orbital ordering and 
the large bandwidth of the Cg states were not considered. It has recently been pointed out by several authors 
that the oxides corresponding to high formal oxidation states may be negative charge transfer systems in the 
Zaanen-Sawatzky- Allen scheme ||^ resulting in an essential modification of the electronic structure, in particular 
in a possible stabilization of an IS states |^ . 

A first physical model to explain the transitions in LaCoOa was proposed and recently revised by Goodenough 
in Q. It was suggested that for trivalent cobalt in LaCoOs the crystal field energy Ac/ is only slightly larger 
than the intraatomic (Hund) exchange energy A^x- Thus the ground state of the Co ion is LS ^Ai (S=0) 
and the excited HS state (S=2), according to §, is only 0.08 eV higher in energy. The corresponding 
one-electron configurations are ^29^^ for the LS state and t2g6g for the HS state. The increase in temperature 
leads to the population of the HS states which is reflected in the magnetic susceptibility measurements. The 
semiconductor-metal transition was interpreted in a rather vague way as a formation of a a* band from the 
localized ionic Cg states. Long-range order of the LS and HS Co ions was assumed based on the results of the 
X-ray diffraction measurements. The recent neutron diffraction experiments |^ did not confirm it, but did not 
exclude a short-range order. 

The theoretical description of LaCoOs is a difficult task, because this system exhibits transition from the 
localized to the delocalized behavior, and the two main existing methods: one-electron band structure calculations 
based on the Local Density Approximation (LDA) and "model hamiltonian configuration interaction" approach 
(which were both used for the compound under consideration 0]) were constructed for the, respectively, 
completely itinerant and fully localized cases. 

Recently a paper of Sarma et al was published [|| where the results of the LSDA calculations for LaAfOs 
(Af=Mn, Fe, Co, Ni) were presented. It was claimed there that LSDA results gave good agreement with the 
X-ray photoemission spectra. Such an agreement is not necessary a proof of the adequate description of the 
electronic structure. A well known example is NiO, where LSDA gives a sharp peak at the top of the valence 
band in agreement with XPS, but the nature of this peak is Ni3d t2g minority-spin character, while it is generally 
accepted that this peak is derived from d?L final state. The real proof of the correctness of the one-electron 
approximation would be an absence of the satellites in the photoemission spectra. Also the XPS spectra, with 
which a comparison was made in were measured at the room temperature, where LaCoOs is already not in 
the low-spin state. So XPS must be compared with LSDA result corresponding to the magnetic solution and not 
a non-magnetic one as in |8| . 

Anisimov et al proposed a so-called LDA+U method |^ which combines in one calculation scheme LDA and 
Hubbard model approaches and is able to treat on a " first- principle" basis the systems with strong Coulomb 
correlation. It was demonstrated that this method (in contrast to standard LDA) could describe the existence 
of the different states of the system under consideration which are close in total energy (holes in doped copper 
oxides ]lO{ and transition metal impurities in insulators [pT[). 

In this paper we present a study of an electronic structure of LaCoOa in the LDA-f U approach. In contrast 
to the standard LDA there are several stable solutions corresponding to different local minima of the LDA-I-U 
functional. We have found a nonmagnetic insulating ground state in agreement with the experiment. We have also 
found two orbitally-polarized magnetic solutions corresponding to IS states (one of them is gapless semiconductor 
and the other is metal) and semiconducting magnetic solution corresponding to a HS state of the system, which 
lies much higher in energy. In addition we found that the Cg states form broad bands whereas the t2g states 
exhibit narrow bands, split by Coulomb interactions into the lower (occupied) and upper (unoccupied) Hubbard 
bands. Using the results of these calculations, we propose an interpretation of the behavior of LaCoOa. According 
to our scheme, there first occurs with increasing temperature a transition from a LS (nonmagnetic) insulating 
ground state to a state with an IS (configuration t^gC^g). Due to strong Jahn- Teller nature of this configuration, 
this state may develop orbital ordering. The orbitally-ordered state turns out to be nonmetallic (actually nearly 
zero-gap semiconductor) in our calculations. With the further increase of the temperature the orbital ordering 
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may be gradually destroyed which can explain the transition to a metallic state observed in LaCoOs at 400 - 
600 K. We hope that this new information will eventually lead to a better understanding of LaCoOa and other 
similar materials. 



2 CALCULATION METHOD 

The LDA+U method Q is based on the assumption that it is possible to separate all electrons in the system 
in two subsets: the localized states (for LaCoOs these are the Co 3d-orbitals) , for which Coulomb intrashell 
interactions are described by the Hubbard-like term, and the itinerant states, where the averaged LDA energy 
and potentials are good approximations. 

It is known that LDA calculations can provide all the necessary model parameters (such as the Coulomb 
parameter U the exchange J, hopping parameters describing hybridization etc.) on a first- principle basis, 
but the one-electron structure of the LDA equations with the orbital-independent potential does not allow to use 
these parameters in the full variational space. LDA-I-U overcomes this problem by using the framework of the 
degenerate Anderson model in the mean-field approximation. In this method, the trial function is chosen as a 
single Slater determinant, so it is still a one-electron method, but the potential becomes orbital-dependent, and 
that allows one to reproduce the main features of strongly correlated systems: the splitting of the d-band to the 
occupied lower Hubbard band and unoccupied upper Hubbard band. 

The main idea of our LDA+U method is that the LDA gives a good approximation for the average Coulomb 
energy of d-d interactions Eav as a function of the total number of d-electrons N — X^ma "mo-, where rima is the 
occupancy of a particular dmCT-orbital: 

But LDA does not properly describe the full Coulomb and exchange interactions between d-electrons in the 
same c?-shell. So we suggested to subtract Eav from the LDA total energy functional and to add orbital- and 
spin-dependent contributions to obtain the exact (in the mean-field approximation) formula: 

E Ej^i^j^ (U ^ tj "^^2 ^ ^ ^mm' ^m(T^m' — (T ^ ^ ^ (Unim' Jvava'^^raa^va' a • (2) 

ra^m' ,CT m^m' ,m' .a 

The derivative of Eq.(2) over orbital occupancy rima gives the expression for the orbital-dependent one-electron 
potential: 

Vmair) = VLDA{r) + ^(J7mm' " Uef f)njn' -cr + ^ {Umrn' " Jmrn' ~ Ueff)nmcr + ^eff{-^ - nma) " -^J ■ (3) 

The Coulomb and exchange matrices Umm' and Jmm' are: 

k 

Jmm' = Y.^kF\ (5) 
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Table 1: Basis set and MT-sphere radii used in the calculation. 



Atom 


kI = -0.01 Ry 


kI = -1.0 Ry 


kI = -2.3 Ry 


Rmt, a 


La 


6s6p5d4f 


6s6p5(i4/ 


6s6p5(i 


1.77 


Co 




AsApM 


AsAp 


1.26 





2s2p 


2s2p 


2s 


0.66 



where the F'' are Slater integrals and {bn\Ykq\lm') are integrals over products of three spherical harmonics 

For d-electrons, one needs F'^,F^ and F* and these can be linked to the parameters U (direct Coulomb 
interaction) and J (intraatomic exchange) obtained from the LSDA-supercell procedures via U — F'^ and 
J = {F^ + F^)/14, while the ratio F"^ / F^ is to a good accuracy constant ^ 0.625 for 2id elements |l^. For 
LaCoOs, the Coulomb parameter U was found to be 7.8 eV and the exchange parameter J=0.92 eV. 

The LDA+U approximation was applied to the full potential linearized muffin-tin orbitals (FP-LMTO) cal- 
culation scheme iQ. Crystallographic data being used in the calculations were taken from ||]. According to 
them, LaCoOs has a pseudocubic perovskite structure with a rhombohedral distortion along the (111) direction. 
The unit cell contains two formula units. Since this rhombohedral distortion is small (the largest rhombohedral 
angle is 60.990° at 4 K), we use the conception of t2g and eg orbitals, as referred to in the cubic setting in the 
following discussion. Temperature was introduced in our calculations only by the change of lattice parameter 
and atomic positions according to the data of Ref.||]. The most detailed description of the technical aspects 
of FP LMTO calculations for the perovskite-type complex oxides can be found in The optimal choice of 

the basis set for describing the valence band and the bottom of conduction band of the LaCoOs is presented 
m Table Since U-correction is applied for the Co d-orbitals, the value of the muffin-tin (MT) radius for Co 
was chosen close to its value in metallic Co in order to get the full 3(i-density inside the sphere. For the correct 
description of the wave functions in the interstitial region, we expanded the spherical harmonics up to the value 
of lmaa;=5, 4, 3 for La, Co and O MT-spheres correspondingly. The Brillouin zone (BZ) integration in the course 
of the self-consistency iterations was performed over a mesh of 65 k-points in the irreducible part of the BZ. 
Densities of states (DOS) were calculated by the tetrahedron method with 729 k-points in the whole BZ. 

As the potentials for the various d-orbitals of Co are different in LDA-t-U, it is not obvious a priori what 
will be the final symmetry of the system under consideration. In preliminary calculations we assumed for the 
simplicity that LaCoOa has a perfect perovskite-type cubic lattice. To allow the system to choose the appropriate 
symmetry by itself, we perform an integration not over 1/48 part of BZ as for cubic symmetry but over 1/8 {T)2h 
symmetry group). The resulting symmetries were found to be cubic Oh for the ground (LS) state and tetragonal 
T)ih for the excited states. Practically the same situation occurs in real R3c symmetry: the occupancies of 
xy, yz, zx orbitals are about the same for the low-spin configuration as well as that of 3^;^ — r^, ~ orbitals. 
The degeneracy of xy and yz, zx orbitals is broken in the other spin states. 

3 RESULTS AND DISCUSSION 

3.1 Homogeneous Solutions 

We start by considering the results of the calculations for homogeneous regimes, without extra superstructure. 
The possibility to get a solution with an orbital ordering is considered in section [3.2[ 

The detailed results of the calculations with crystallographic data corresponding to 4 K are presented in Table 
||. One must start with the equal occupancies of all three t2g orbitals and also of two eg orbitals (e.g. i^gCg) 
to obtain the LS configuration. The initial spin polarization vanished during self-consistency iterations resulting 
in a nonmagnetic final solution. The charge excitation spectrum has a semiconducting character in accordance 
with experimental data ||l^ (and in contrast to the LDA result 0) with the energy gap equal to 2.06 eV. The 
top of the valence band (Fig.l) is formed by the mixture of oxygen 2p states with Co t2g orbitals and the bottom 
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Table 2: Total energy difference, character of energy spectrum, occupancies of various d-orbitals of Co and 
magnetic moments for various Co spin-states in LaCoOs. 
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yz 
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total 








2.06 


T,i 


0.98 


0.98 


0.99 


0.32 


0.33 


7.20 





low-spin 
























0.24 


0'' 
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0.98 


0.98 


0.98 


0.84 


0.84 


7.13 


2.11 


intermediate-spin 






i 


0.07 


0.98 


0.99 


0.23 


0.24 








0.65 


0.13 


T 


0.99 


0.99 


0.99 


1.00 


1.00 


6.78 


3.16 


high-spin 






i 


0.99 


0.08 


0.11 


0.33 


0.30 







" Total energy difference relative to the energy of iig^g configuration. 
Half-metal. Total magnetic moment per unit cell is 2 fis- 



of conduction band predominantly by the Cg orbitals. 

From the X-ray photoemission and absorption spectra 0| the gap was estimated to be 0.9 ± 0.3 eV. The 
optical measurements gave w 0.5 eV. The larger value of the calculated gap can be explained by the well 
known fact that a mean-field approximation, which is the basis of our LDA-I-U approach, usually overestimates 
the tendency to localization and hence the values of the gap. 

The aim of our work was to find solutions corresponding to the higher-spin states of Co ions. Usually the HS 
state is described as a ilg^g configuration with the maximum value of spin S=2 (magnetic moment ^(70=4 /i^). 
This corresponds to the purely ionic model, and hybridization of the Co 3d-orbitals with the oxygen 2j3-orbitals 
and band formation in a solid can renormalize significantly this ionic value. Such kind of renormalization was 
obtained in our calculations (see Table The initial f^g^^ configuration (two holes on d^^, d^^ orbitals of the 

t2g set and two on d^^2_j.2, d^2_y2 of the Cg set with spin-down (minority) spin projection) results in the self- 
consistent solution with a magnetic moment of fj,d^co—3.l() fJ-B- The total energy of this HS solution is 0.65 eV 
per formula unit higher than the ground state LS ^fo^g configuration. The unexpected result is that there exists 
an intermediate-spin solution (second line in Table magnetic moment value /^c;-Co=2.11 /is) which is lower 
in total energy than the HS solution (0.24 eV per formula unit relative to the LS ground state). This solution 
was obtained when we started from the initial configuration ^Ig^gi where only one electron was transferred from 
the t2g to Cg states. The final self-consistent configuration is better described as a state with a partially 

filed CT*-band (i^gf*) with the occupancies of the d'i^^2_^-2 and d^^2_y2 orbitals equal to 0.84. In a configuration 
interaction language used in the cluster calculations this may be compared to a d^+d'^L state where L denotes 
a hole on the oxygen. 

The IS state in our calculations turns out to be metallic (see, however, section 3.2), and HS state is semi- 
conducting (Fig. 2, 3). The reason for this is that the antibonding a* band (formed by eg orbitals) is very broad 
and the band splitting is not strong enough to create a gap in the case of IS state. When the band becomes 
completely filled (HS state), a small gap between e| and e^ bands appears. 

In Table ^ the occupation of different orbitals in different configurations is presented. One notices that the 
actual occupation is different from the formal " chemical" one and corresponds on the average not to 6 but rather 
to 7 electrons in d-shell. The reason for this is the strong hybridization of the empty eg orbitals with oxygen 
2p-orbitals. In the LS ground state every eg orbital, which is formally empty, has actually an occupancy of about 
0.33 resulting in 1.3 additional d-electrons above the formal six-electron configuration. In the case of excited IS 
state we have partially filled broad (T*-band with the total number of Cg electrons increased by 0.85 as compared 
to the LS state and the number of t2g-electrons 0.92 less than in LS state with as a result practically unchanged 
total number of d-electrons. Therefore, we used the notation t2gCr* in the Table ^ to prevent misunderstanding 
and in accordance with notation proposed in ^ . One may also say that despite the formal oxidation state Co^+ 
the real configuration e.g. in the IS state is a mixture of i^g^g and t2ge^L configurations. 

Some other points concerning the results obtained are worth mentioning. As one can see from Fig.l, the 
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Figure 1: The total and partial densities of states obtained in LDA+U calculations for LaCoOa with Co ions in 
low-spin state (^23^° configuration). For d-Co partial density of states solid line denotes t2g orbitals and dashed 



line - Cg orbital. 
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Figure 2: The total and partial densities of states obtained in LDA+U calculations for LaCoOs with Co ions in 

intcrmcdiatc-spin state {t2g<^* configuration). Fermi level is denoted by vertical dashed line. Arrows corresponds 
to spin- up and spin-down spin projection. The other denotations are the same as on Fig.l. 
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Figure 3: The total and partial densities of states obtained in LDA+U calculations for LaCoOa with Co ions in 
high-spin state {tigSg configuration). Denotations are the same as on Fig.l and Fig. 2. 
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photoemission holes in a LS state are mainly formed by oxygen states (note that in Figs. 1-3 the partial densities 
of states of constituent atoms are given per one atom, and there are 3 of oxygen atoms per formula unit). 
Consequently the final state of the XPS for LS state formally corresponds to d^L configuration. At the same 
time the electronic excitations mostly reside on Co Cg orbitals hybridized with p-states of oxygen. However, in 
the IS case which in our calculations is metallic (Fig. 2), the states at the Fermi level contain comparable weight 
of both Co Eg and O p orbitals. 

The IS state obtained is fully polarized which is due to the half-metallic ferromagnetic nature of the solution: 
the magnetic moment per formula unit is 2 i-e. it corresponds to the spin S=l. In that sense the calculation 
agrees with the ionic picture, in which the IS configuration of Co is d^ (^ig^g); with S=l. 

The IS state of LaCoOa in our results turned out to be half-metal ferromagnet, in which only the electrons 
of one spin projection are at the Fermi level (the other spin subbands are fully occupied). Ferromagnetic arrange- 
ment of the magnetic moments may be of course an artifact of the scheme used. In principle, LDA-I-U method 
is intended to describe the local correlation effects such as formation of the Hubbard gaps and local magnetic 
moments, which should still be present in the paramagnetic phase. However technically in the calculations we 
have to assume certain long-range magnetic order in accordance with translation symmetry of the crystal. What 
the results obtained in such a way do really tell you, is whether there is a gap or not and what is the value of 
the local magnetic moment. 

Although the IS state in the calculations turns out to be metallic, the density of states at the Fermi level 
is very low, n(Ei?)=0.36 states per eV per one formula unit. This probabl y ind icates that it would not be so 
difficult to make this state semiconducting, see the discussion below (section 3.2). 

The results presented above show that the first excited configuration - that with the IS - lies only 0.24 eV 
higher than the ground state with the LS. Experimentally it is established that there is a transition from LS 
nonmagnetic state to the magnetic one with the increase of temperature. According to the Ref.[|| this transition 
occurs in the vicinity of 90 K. As the closest magnetic state is that with the IS (HS state lies, according to 
our results, much higher at 0.65 eV), we ascribe nonmagnetic - magnetic transition in LaCoOs to the LS - IS 
transition. 

It is well known (see, e.g., [Q) that Co in the HS state has larger ionic radii than in LS one, and LS - HS 
transitions are accompanied by the increase of the volume (and vice versa). Keeping that in mind we carried 
out the calculations of the electronic structure of LaCoOa at different lattice parameters which can imitate the 
influence of the temperature (via thermal expansion). Fig. 4 demonstrates the values of total energies for various 
spin states of Co relative to the energy of ^ig^g state at the lattice parameter corresponding to 4 K. In this figure 
actual lattice parameters used in calculations are shown on the horizontal axis together with the temperature 
scale (we used the lattice parameters as a function of the temperature measured in Ref. ) . One sees that with 
the increase of the specific volume (or with the increase of temperature) the energy of the IS state crosses that of 
the LS state, which corresponds to the LS - IS transition. In our calculation this crossover occurs at about 150 K 
- somewhat larger that the experimental value about 90 K. Nevertheless this transition temperature is correct 
by order of magnitude which for such ah initio calculations is rather satisfactory. Note that the HS state still lies 
high enough even at the temperatures about 600 K although we can not exclude that it could have become the 
ground state at still larger specific volumes or temperatures. 

The total energy for all three solutions (LS, IS, and HS) decreases with the increasing volume and minima 
are achieved only for lattice parameters corresponding to high temperature. At this value of the volume the 
total energy of the IS state is lightly lower than LS. It is well known that the equilibrium values of volume 
calculated with LSDA are always a few percent less than the experimental volume, because LSDA overestimates 
cohesion. LSDA-fU, on the other hand, underestimates cohesion due to its treating the d states as localized. 
In both cases the calculated lattice parameters significantly deviate from the experimental values. Goodenough 
suggested that the magnetic transition in LaCoOa is caused by the fact that the crystal-field energy 5cf is only 
slightly larger than the intra-atomic (Hund) exchange energy 5eX. The crystal-field energy is determined by the 
3d-2p hopping parameters, which strongly depend on crystal volume. Hence if one will perform calculations with 
lattice parameters corresponding to the equilibrium volume, then the delicate balance between the Scf and SeX. 
Our results show that for low-temperature lattice parameters is larger Scf than the S^x, but already a small 
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Figure 4: The total energies for various spin states of LaCoOs relative to the energy of ffg^g state at 4 K versus 
R3c lattice constant. The correspondent temperatures are marked also. 
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Figure 5: Scheme representation of various Co d^+dJL configurations in low- (a), intermediate- (b) and high- 
(c) spin states. Open circle denotes a hole in oxygen p-shell. 

increase on volume (corresponding to increasing the temperature to 150 K) would reverse this ratio. 

Thus from the results obtained we conclude that the nonmagnetic - magnetic transition in LaCoOs most 
probably occurs not between LS and HS states, but between low- and intermediate-spin states. The reason for 
the stabilization of the IS state may be understood from the following consideration. Different configurations of 
Co'^+ {d^) are illustrated in Fig. 5. Here and Cg are atomic c?-states splitted by the crystal field (the splitting 
A=10 Dq). In purely ionic picture one should expect that, depending on the ratio between A and intraatomic 
exchange interaction J either LS or HS states should be stable (if the energy of the LS state Els is taken to be 
zero, then E/5=A-J and E//5=2A-6J, so that LS state is the ground state for A >3J and HS state - if A <3J; 
IS state in this scheme always would lie higher). However, as we mentioned in the Introduction, the oxides with 
unusually high valence of transition metals have a tendency to have the d-shell occupancy corresponding to a 
lower valence, with the extra hole being predominantly located on oxygen (one gets a configuration d'^L instead 
of d^). In this situation the d-p hybridization is especially strong and plays a crucial role. These considerations 
are also supported by our calculations, see Table 

Especially important is the hybridization with the Cg orbitals. One can show that it is stronger for the IS state 
than for the HS one: there are more channels of p-Cg hybridization in this case (the most favorable configurations 
of the type d'^L, admixed to d^, are also shown in Fig. 5). In other words, if these configurations (d^L) would be 
dominant, then the state corresponding to the ground state of the configuration would have lowest energy - 
and it is known that for d'' (Co^+) it is just the configuration of Fig. 5b - i.e., the one giving a total state of IS. 
That is the physical reason why the IS configuration lies below the HS one and may become the ground state for 
an expanded lattice (higher temperatures). 

The nonmagnetic - magnetic transition occurring at ^^90 K which we now associate with the LS - IS transition, 
is described quite reasonably by our model. The transition temperature obtained in our calculations (~150 K) 
lies not so far from that experimentally observed. In Ref. |^ the relative change of the lattice parameter ( Aa) at 
this transition was fitted by the relation (see ||Tg||) 

Aa = " *^^Pi-Eq/'^'^) eg-) 
l + iy*exp{-Eq/kT)' ^ ' 

where Eg is the energy gap (it was taken in |^ as a fitting parameter), and v is degeneracy of the magnetic 
state. As the high-temperature magnetic ground state was believed to be the HS state in Ref. ||] , the degeneracy 
was taken v=\b (triply degenerated t2g orbital times the spin degeneracy (2S-I-1) with S=2, see Fig. 5c). In our 
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new interpretation the degeneracy will be iy=18 (triply degenerated t2g orbital times double degeneracy of Cg 
orbital times (2S+1) with S=l, see Fig. 5b). Thus the fit of Aa would be as successful with this interpretation 
as with the one given in ||^, with actually nearly the same value of the energy gap Eg. 

The picture of an IS state can also help to resolve one more problem mentioned in Ref. |2| : that the correlations 
between magnetic sites in LaCoOa are not antiferromagnetic but (weakly) ferromagnetic. One should expect only 
antiferromagnetic correlations between high-spin Co^~^ ions in which the e^-shell is half-filled (e^). However in 
IS Co ions have nominally configuration, and especially if these Cg orbitals are ordered (see discussion below) 
one can have a ferromagnetic exchange interaction (the well-known example is e.g. ferromagnetic K2CUF4 pO[|). 

3.2 Possible Orbital Ordering in LaCoOs 

The conception of LS - IS transition describes quite reasonably the nonmagnetic - magnetic transition observed 
in LaCoOa. At the same time it was established that the first transition at ~90 K leaves this compound 
semiconducting whereas our calculated IS has a metallic character of the energy spectrum. Is it possible to 
overcome this disagreement? 

As is clear from Fig. 5b, the IS state has a strong double-degeneracy (e^ configuration). For localized electrons 
it would be the typical Jahn- Teller situation. In particular for KCuFs it was proposed in Ref.pO|] and it was 
confirmed by calculations in Ref. pT[ | that there exist the special kind of orbital and magnetic ordering of d-ions in 
simple cubic lattice with one electron or hole on doubly-degenerate eg level. The situation in IS LaCoOs is very 
similar to the case of KCuFs (keeping in mind that there exist permanent rhombohedral distortion and there are 
already two Co ions in an elementary cell in LaCoOa, the "antiferro" orbital ordering may be consistent with 
the existing structural data). 

We checked this possibility by repeating the calculations, assuming now that there exist an orbital ordering. 
The structure of the corresponding ordered IS state was taken as consisting of ferromagnetic planes (001) whereas 
the direction of spins is opposite in the nearest planes. The occupied Cg orbitals were assumed to alternate; as 
a starting point the orbitals in two sublattices were taken as drj.2_z2 and dy-2_^2, see Fig. 6. To supplement this 
Figure with the information about the other d-orbitals, let's consider the Co ion in t^gcj configuration which has 

d}y2_^2 occupied eg orbital (corresponding site is marked as 1 on Fig. 6). For this state the other three eg orbitals 
are empty and all the t2g orbitals are occupied, with the exclusion of d^^ o^'^. This choice minimizes the Coulomb 
interaction energy of t2g and eg electrons. For the neighboring Co ion (site 2 in Fig. 6) the eg electron is placed 
in d\_ 2 orbital and t2g hole is placed in d^^ orbital. Calculations were made for the real rhombohedral crystal 
structure with lattice parameters corresponding to 71 K. With this state as a starting point the self-consistent 
calculation was now carried out. The resulting densities of states corresponding to such orbital ordered IS state 
are presented in Fig. 7. 

First of all it is necessary to point out that such kind of spin and orbital ordering is stable with respect to the 
self-consistency iterations and hence there exists the solution with such symmetry. The character of an energy 
spectrum turned out to be gapless semiconducting-like. The wide eg band which crossed the Fermi level in the 
case of IS state without orbital order (see Fig. 2) is now splitted. The top of the valence band is formed (for one 
sublattice, e.g. for Co site 1) by dy2_^2 states and the bottom of conduction band - by d^x^-^^ states (Fig. 7b, c) 
hybridized with the oxygen states. As to the t2g densities of states, the character of their energy distribution 
is practically not changed in comparison with the case of IS state without orbital ordering. The value of spin 
magnetic moment of Co is 1.87 /is in this orbitally ordered IS state and it is less than for one without orbital 
order (2.11 hb)- 

When we compare the two IS states with the fixed spin magnetic moment 1.87 fiB, the total energy of the 
orbitally ordered IS state is 0.11 eV lower. The necessity of using a fixed magnetic moment is due to the limitation 
of the band structure scheme. The magnetic state of LaCoOa is known to be paramagnetic (with the possible 
short-range spin or orbital ordering). At the same time as we already mentioned above, it is necessary to set 
a long-range magnetic order in calculations (antiferromagnetic for orbitally ordered IS state and ferromagnetic 
for one without orbital order). This long-range magnetic ordering imposed (in addition to the intraatomic 
exchange) influences somehow both the local magnetic moment and the total energy. In the case of long-range 
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Figure 7: The total and partial densities of states obtained in LDA+U calculations for LaCoOs with Co ions in 
ordered intermediate-spin state, a: Total density of states (per 4 formula units and both spins); b: partial Co 
eg densities of states e.g. for Co site 1 in Fig. 6: solid line - for (iy2_j,2 states, dashed line - for d^x^-r^ states; c: 
the same as in b in the vicinity of the Fermi level. Denotations are the same as on Fig. 2. 
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antiferromagnetic order the value of the magnetic moment is apparently underestimated in comparison with the 
real paramagnetic state, for the ferromagnetic one - overestimated. Thus to estimate the relative stability of 
different phases we compare the energies calculated for the same value of the magnetic moment. Our calculations 
show that in the case of equal spin values the orbitally ordered IS state is more preferable than the one without 
orbital order. 

One can interpret now the first nonmagnetic - magnetic transition (at about 90 K) as a transition of Co 
ions from LS to IS state with the specific orbital ordering of occupied eg orbitals. Our results show that this 
transition occurs at temperatures lower than 150 K. Under this transition a magnetic moment on Co sites appears 
whereas the material is still nonmetallic. The spin-only value of the effective magnetic moment, in the case of IS 
state (S=l), is /Lte//=2.83 /is - somewhat below experimental value of 3.1-3.4 fiB- It will be increased by the 
orbital contribution, because there is in principle an unquenched moment of t2g electrons for the configuration 
^2g^l (which e.g. for Co^+ ions the increases the value of /ie// typically by '--^0.3 /is). We cannot treat this 
effect numerically because our codes do not account for the spin-orbit interaction; however we expect on physical 
grounds that this effect should be present (which, by the way, would change somewhat the occupation of t2g 
orbitals). An extra change of fieff may be due to the very strong d-p hybridization which leads to a large 
contribution to the total wave function of a state dJL (ilgCgL). In this state the effective moment of Co ion is 
that of the high-spin Co^+ (S=3/2) compensated by the opposite polarization of the oxygen p-shell (S=-l/2) 
(note that this picture is supported by the calculations (see Table |^) where local Co spin moment is 2.11 /is and 
the total spin moment per unit cell is 2 fis)- This fact may significantly change the value of fJ-eff measured in 
the susceptibility experiments. 

The IS state with the nominal configuration ^23^9 may have an orbital ordering because of the doubly de- 
generate eg orbitals (strong Jahn- Teller nature of this configuration). This turned out to be the case in our 
calculations: the orbitally ordered state may be crudely described as an alternation of the occupied (1^2 and 
(1^2 _y2 orbitals. 

According to our calculations, this ordering leads to the splitting of the eg ((7*)-band, leading practically to 
the zero-gap situation, with the Fermi level lying in the gap. The second gradual transition to the metallic state 
('^GOO K) with the increase of the effective magnetic moment value up to fieff—'^-O fiB may then be associated 
with the disappearance of the orbital ordering with increased temperature still within the IS state. This IS state 
without orbital order is calculated to be a metal with a larger magnetic moment. Our calculations show that the 
second transition has no relations to the HS state of Co ions because the HS solution lies still high in energy and 
if realized would have been semiconductor. 

4 CONCLUSION 

The calculations of the electronic structure in LDA+U approximation for LaCoOa were performed. At 4 K 
the lowest total energy has the nonmagnetic insulating solution with Co ions in LS state. Three excited state 
configurations (two IS and one HS) are also found to be stable, with local magnetic moments on Co sites. The 
t2g states exhibit narrow bands whilst the eg states form broad bands. With the increase of lattice parameter 
corresponding to the thermal lattice expansion two transitions can be expected: the first occurs from LS state 
to IS state with the orbital ordering which in our calculation is a zero-gap semiconductor. The second transition 
occurs within the IS state and is connected with gradual disorder of occupied eg orbitals. We believe that 
although many details are not clear yet, the concept of the nonmagnetic - magnetic and semiconductor - metal 
transitions in LaCoOa as connected mainly with IS states is a good candidate to explain the properties of this 
interesting material which remains a puzzle for so long. 
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